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Abstract 

The second-order Mpller-Plesset ab initio electronic structure method is used to compute points 
for the anharmonic mode -coupled potential energy surface of V-methylacetamide (NMA) in the 
trans c t configuration, including all degrees of freedom. The vibrational states and the spectroscopy 
are directly computed from this potential surface using the Correlation Corrected Vibrational Self- 
Consistent Field (CC-VSCF) method. The results are compared with CC-VSCF calculations using 
both the standard and improved empirical Amber-like force fields and available low temperature 
experimental matrix data. Analysis of our calculated spectroscopic results show that: (1). The 
excellent agreement between the ab initio CC-VSCF calculated frequencies and the experimental 
data suggest that the computed anharmonic potentials for /V-methylacetamide are of a very high 
quality. (2). For most transitions, the vibrational frequencies obtained from the ab initio CC- 
VSCF method are superior to those obtained using the empirical CC-VSCF methods, when com- 
pared with experimental data. However, the improved empirical force field yields better agreement 
with the experimental frequencies as compared with a standard AMBER-type force field. (3) The 
empirical force field in particular overestimates anharmonic couplings for the amide II mode, the 
methyl asymmetric bending modes, the out-of-plane methyl bending modes, and the methyl distor- 
tions. (4) Disagreement between the ab initio and empirical anharmonic couplings is greater than 
the disagreement between the frequencies, and thus the anharmonic part of the empirical potential 
seems to be less accurate than the harmonic contribution. (5) Both the empirical and ab initio 
CC-VSCF calculations predict a negligible anharmonic coupling between the amide I and other 
internal modes. The implication of this is that the intramolecular energy flow between the amide 
I and the other internal modes may be smaller than anticipated. These results may have important 
implications for the anharmonic force fields of peptides, for which V-methylacetamide is a model. 
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I. Introduction 


There is long standing interest in the spectroscopy and dynamics of Af-methylacetamide (NMA). 
This is due in part because NMA can serve as a simple model for the amide bond in peptides and 
proteins. The importance of the amide group lies in it’s contribution to both intramolecular back- 
bone hydrogen bonding and solvent-protein hydrogen bonding. An experimental estimate for the 
energetics of main chain hydrogen bonding is between 1.3 +/ — 0.6 and 1.5 +/ — 0.3 kcal/mol. 1,2 
In order to accurately represent the physics of the amide bond, it is important to first understand 
both the intramolecular forces and the coupling this amide group undergoes. This implies a detailed 
investigation of the molecular forced fields used for Af-methylacetamide, from both an empirical 
and an ab initio perspective. However, any such study should not only include the energetics of 
the different conformers of NMA, but also anharmonic vibrational spectroscopy and the effects of 
intramolecular vibrational mode coupling. One interesting example of amide coupling may be il- 
lustrated in the two broad subbands found in the vibrational Raman spectrum for the amide I mode 
of NMA in H 2 0? In this work, Chen et. a/ found that these two subbands were probably not 
due to different conformations of the NMA, but could in fact be attributed to coupling between the 
amide I mode to the bending mode of the solvent H 2 0 molecules. However, in the femtosecond IR 
pump-probe spectra of NMAD (Deutrated NMA) Hamm and co-workers found the same doublet 
structure in the amide I band for NMAD in D 2 0 and suggested that this subband could not be due 
to coupling to the OD bending mode of D 2 0. A Along a similar line of work, Pajcini and Asher 
investigated excitonically coupled states which mix the primary and secondary amide transition 
of N-acetylglycinamide (NAGA). 5 By looking at the polarized-oriented Raman spectra of single 
crystal NAGA, they were able to show that the primary and secondary amide III and amide I vibra- 
tions are excitonically coupled. Therefore it is still an open question if the coupling of the amide 
I mode is to other internal modes or to the solvent modes. The question of intra/intermolecular 
mode coupling becomes even more important when considering vibrational energy redistribution 
pathways as measured by femtosecond nonlinear-infrared spectroscopy. 4 

Notwithstanding the deep and detailed experimental 6-13,3,14,6 and theoretical 15 22 investigations 
of N-methylacetamide, relatively little is still known about the intramolecular coupling between 
different vibrational modes and the effect of this on spectroscopy. In previous ab initio theoreti- 
cal calculations of the energetics and harmonic vibrational spectra of N-methylacetamide, it was 
illustrated that one needs to scale the harmonic frequencies to reasonably compare with the ex- 
perimental data. 16 In addition, the different conformations of NMA (cis and trans ) require slightly 
different scaling factors. Furthermore, when the amide group is hydrogen bonded to water, this also 
requires rescaling of the force constants. 16 Thus, although the ab initio theoretical calculations of 
N-methylacetamide were of a reasonable level, the calculated harmonic vibrational frequencies re- 
quired scaling in order to compare well with known experimental data. 16,22,17,21,18 Because scaling 
is an empirical and non-transferable procedure, vital information regarding the spectroscopy and 
more importantly the anharmonic nature of the underlying molecular force field, is lost. Therefore, 
it is beneficial to derive an ab initio molecular force field that will not only directly reproduce 
the vibrational spectroscopy but also the anharmonic intramolecular mode coupling nature of the 
molecule. 

One can, in principle, measure vibrational frequencies quite accurately with low tempera- 
ture vibrational spectroscopy, either in rare gas matrices, 23 in a nozzle jet expansion, 24-29 or in 


2 



He droplets. 30-32 This is because at lower temperatures the effects of line broadening and fre- 
quency shifting will be minimized. For example, in the recent superfluid He droplet experiments 
of Huisken et. al, the authors could spectroscopically resolve different conformers of glycine 
based on the OH stretching frequency. When the OH in glycine is internally hydrogen bonded, 
the frequency will be red shifted by 300 cm -1 . 23 Furthermore, recent developments in the field of 
two-dimensional infrared spectroscopy show substantial promise in the ability to not only ascertain 
anharmonic shifts in vibrational frequencies, but also to resolve mode -coupling effects. 33,34,4,35-38 
This newly emerging experimental field requires a companion in the theoretical understanding of 
anharmonic mode-coupled vibrational spectroscopy which we shall begin to address by employing 
the Correlation Corrected Vibrational Self-Consistent Field method. 

When the potential is available in an analytical form, one can perform a Vibrational Self- 
Consistent Field (VSCF) calculation in order to determine the frequencies and relative degree of 
mode coupling in the system. This has proven to be quite accurate and useful for a wide variety of 
problems such as (. 4 r) 13 clusters, 39 peptide-water complexes, 40 glucose, 41 and BPTI with almost 
200 hydration waters 42 Moreover, by extending the VSCF method using perturbation-corrections 
treatments, Jung and Gerber were able to deal quite successfully with the highly anharmonic cou- 
pled mode system, ( H 2 0 ) n , n < 8. 39 Obviously, this method is as successful as the underlying 
analytical force field is accurate. Until recently, only empirical force fields were used in the calcu- 
lation of spectroscopy at an anharmonic level for all but very small molecules. Recently, Chaban 
et al. have developed an algorithm that combines anharmonic vibrational spectroscopy with direct 
calculation of ab initio potentials. 43 The method accounts for the anharmonicities and couplings 
between vibrational modes using the Correlation Corrected Vibrational Self-Consistent Field (CC- 
VSCF) approach and, therefore, provides a superior alternative to the existing techniques based 
on scaling the harmonic vibrational frequencies. This method does not require a fitting of an an- 
alytic potential function, nor high order derivatives. Thus, direct calculation of the anharmonic 
vibrational spectroscopy is feasible for molecules up to 15 atoms. The direct ab initio VSCF 
method was successfully applied to the calculation of the fundamental excitations of a number of 
hydrogen -bonded systems, 44 different conformers of glycine, 40 and a glycine-water complex. 46 It 
is our aim to apply this method for the calculation of the anharmonic vibrational spectrum of trans 
N-methylacetamide. In particular, we shall investigate the intramolecular mode coupling between 
different pairs of vibrational modes, as compared to the results obtained from a calculation of the 
anharmonic vibrational spectrum using an empirical force field. We shall be able to directly as- 
sign a coupling strength to each pair of modes. This work represents a continuing investigation 
into the nature of anharmonic coupling in biological molecules, and can be directly compared to 
experimentally determined two-dimensional vibrational spectrum. 

The paper is organized as follows: A general discussion of the vibrational self consistent field 
method is given in section III. Section III, part A will outline the empirical force field and the 
empirical VSCF method, while part B will discuss the ab initio VSCF method and the energy 
calculations. The results and a discussion is given in section IV and the conclusions follow in 
section V. 


II. Vibrational Self Consistent Field Method (VSCF) and Correlation-Corrected 
VSCF 

Consider the vibrational Schrodinger equation: 

{-% £ + v (Qi. •••> <2<v)]*M., Qs) = £*(<?. Q») 

l jz=l oQj 

2.1 

treated here in mass-weighted normal mode coordinates, where Q } is the j th mass-weighted nor- 
mal mode coordinate and V(Q U ..., Q N ) is the full potential which includes anharmonicity and 
coupling between all of the modes. One may simplify equation 2.1 by utilizing the separable 
Hartree approximation, which is the basis of the VSCF method: 47-50 '» 

^(Qi, Qn) = II ^k(Qk) 

fc=i 

2.2 


This then leads to N single mode VSCF equations of the following form: 

[ -f^f +v '‘ w ‘ )-£kWW * )=0 

2.3 


where the effective potential for mode k is given by: 

v?(Q k ) =< n ^(QdiviQi, Qjv)i n w(Qi) > 

Ijik l^k 

2.4 

Equations 2.3 and 2.4 for the single mode wavefunctions, energies, and effective potentials must 
be solved self-consistently. The VSCF expression for the total energy of vibrational state n is thus 
a sum of all the individual mode energies minus a term which accounts for the double counting of 
the interactions in the energy calculation and has the following form: 

E & U = - (*- 1) < n «iv'(Qi,...,q*)i n«> 

Jt=l k = 1 k = 1 

2.5 

In order to account for correlation effects between modes, an effective perturbation treatment anal- 
ogous to Mqller-Plesset method (MP2) for electronic structure calculations, is employed. In this 
treatment, the expression is up to second order perturbation in AV , and is as follows: 
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k=l 


The second-order correlation-corrected approximation for the energy is: 

F cc _ p vscf . v n; v =1 *;»(&)> P 
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In equation 2.7 we assume that there are no degenerate excited states (n,m), with energy E® — E^. 
The theory behind the correlation-corrected VSCF (CC-VSCF) approximation can be found in 
references 39 and 43-49. 


A. An Expansion Approach to the CC-VSCF Method and the Empirical 
Force Field 


In the present study we utilized an empirical potential which is based on an Amber type of force 
field. Amber is one of the more widely used force fields in many biological studies./citeWeiner84,Weiner86 
The functional form for the analytical potential is as follows: 51 


V n 


V{R) = Y, K r( r - r e q ) 2 + E K »( e - M 2 + E + COS ( n< £ - 7)] + 


bonds 

n Aj 

i<j 


angles 


dihedrals 
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^wr4 )+d '^ 


V 


2.A.1 


Here R is a vector containing the x,y,z coordinates of all the atoms, r represents the norm of the 
bonding vector between any 2 connected atoms, and 9 is the angle between any 3 bonded atoms. 
The dihedral angle, <t>, includes both proper and improper torsion, where 7 is a phase factor and V n 
is the torsional barrier parameter. The last two terms in equation 2. A. 1 represent the Van der Waals 
and electrostatic interactions, which are calculated between atoms of either different molecules, or 
on the same molecule but separated by at least 3 bonds (un = 8 and e/14 = 1-2). Here Rij is the 
distance between atoms i and j, Aij, Bi j are the Van der Waals parameters, and qi is the partial 
charge on atom i. This force field is based on the ab initio work of Cornell et al. and is adjusted to 
fit the experimental vibrational spectrum of trans A/-methylacetamide. 13,16 This was accomplished 
by binning both the torsional parameters as well as the partial charges on the methyl group atoms. 
A Monte Carlo search was performed on the adjustable parameters, in conjuction with a VSCF 
calculation of the vibrational frequencies, in order to reach better agreement with the experimental 
data. The best fit molecular parameters for the implementation of the improved empirical force 
field can be found in Table I. 
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Assuming the system is not extremely anharmonic we can apply a power expansion about an 
equilibrium position to obtain the VSCF effective potential, Vk{Qk) of equation 2.4. Here V k {Qk ) 
is calculated by averaging over all the other l k modes as follows: 

gg QU |g « + ~Ql + 5 E gg Ql < MQd\Qi\MQi) > + 


1 iV n3 Y 

2 E < MQiMMQ:) > +i < UQMMQi) 


> 


2.A.2 

Here, the effective potential is based on an expansion (Taylor Series) up to fourth order about 
a local/global minimum, however a grid-based method is also acceptable. In our present calculation 
we have neglected the cubic and higher order terms with non repeating indices ( QkQiQm , QkQiQmQn ) 
as they will be exceedingly small and in fact vanish in a perturbation treatment which uses the 
harmonic approximation as the unperturbed Hamiltonian. 52 For n > 2 the expansion coefficients, 
d n V 

— — , were calculated using finite differences. 53 The quartic expansion (eq 2.A.2), suffices for such 

oQk 

systems, if one studies low-lying vibrational states, and if the molecule is only moderately anhar- 
monic. However, some systems do exhibit strong anharmonic effects, as illustrated for ( H 2 0) n 
clusters, 54 in which case the expansion method fails completely. For this system, a grid-based 
method is superior to an expansion-based method. In the case of NMA, numerical test indicated 
that the expansion method and the grid method were equivalent, however the grid based method 
requires less cpu time, and thus was chosen. 

In our present calculation, we solved each k VSCF Hamiltonian (eq. 2.3) using the collocation 
method of Yang and Peet. 55 For an adjustable grid, sensitive to each mode, we choose to work in 
the dimensionless variable, qk, where 


Qk - 


Qk 

(A*)< 


2.A.3 


Here A* is the k — th eigenvalue of the Hessian. Convergence was determined when [A^[£*]] < 

k 

0.0001 kcal/mol from one iteration to the next. 


B. Ab Initio Vibrational Self-consistent field Method: 

The calculation of the effective potential, equation 2.4, in an ab initio method will scale unfa- 
vorably (in cpu time) if one were to apply a Taylor Series expansion about a given configuration. 
However numerical test indicated that the pair-interaction approximation used in the ab initio cal- 
culation, and the expansion method (equation 2. A. 2) used for the empirical force fields, were nearly 
equivalent. Therefore, the approach taken was to assume that the potential function can be well 
represented by including interactions between only pairs of modes, and calculating the ab initio 
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potential on a two-dimensional grid. 39,54 Our experience with this approximation, which neglects 
direct coupling between triplets of normal modes, has been very encouraging. 43-45 In this pair-wise 
coupling approximation, the effective potential becomes: 


N 


Vvscf(Q l) •••) Qn ) — ^ V j(Qj ) -l- ^ V,j(Q,, Qj) 

j = i 


2.B.1 


where 


V J (Q J ) = V(Q l =Q,...,Q j ,...Q N = 0) 


and 


Qj) = V(Qi = 0, ...Q u Q v ..., Qn = 0). 


2.B.2 

We stress that our numerical test for several biological molecules have shown that the pairwise 
approximation seems valid for such molecules, and that adding additional mode-couplings does 
not make a significant contribution in terms of the calculated frequencies. 

An electronic structure calculation is first used to calculate the equilibrium configuration. Then 
from diagonalization of the Hessian at this configuration, the normal mode coordinates are ob- 
tained. In order to calculate the potential at each normal mode point, a transformation back to 
Cartesian coordinates is made, and the potential energy is then calculated at this displaced geome- 
try. The diagonal and mode-mode pair-coupling potentials of equation 2.B.2 were calculated on a 
8-point and 8X8 grid which were later interpolated to a 16-point and 16X16 grid. 

Finally, the IR intensities are calculated using the dipole moment estimated along the normal 
coordinate the VSCF wavefunctions for the ground and excited vibrational states as: 44 


U = > I 2 - 

2.B.3 

In equation 2.B.3, uij is the CC-VSCF vibrational frequency for mode i and Qi ), are 

the ground and excited state VSCF wavefunctions. 

All ab initio calculations were performed using the GAMESS electronic structure program, 56 
with second-order Mpller-Plesset perturbation theory (MP2). The basis set used was a Dunning- 
Hay double zeta + polarization (DZP). 57 This level of theory was found to give satisfactory results 
for previous ab initio-WSCF calculations on small biological molecules: glycine, 44 and a glycine- 
water complex. 46 The number of ab initio single-point energy calculations required for both the 
diagonal and pair-coupling potentials is: 
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2.B.4 


where N mode is the number of normal modes, which in this case is 30, and N grid is the number of 
grid points in one dimension. For the present calculation the number of points, on an 8-point grid, 
is 28,080. The total ab initio VSCF calculation required approximately 4 months of CPU time on 
a Octane SGI workstation. 

III. Results and Discussion 

In this work, anharmonic vibrational states and frequencies of W-methylacetamide were ob- 
tained for three potentials; a standard AMBER force field, an improved empirical force field cal- 
culated to agree with existing vibrational frequencies of this system, and the MP2/DZP ab initio 
anharmonic potential energy surface. For each of these potentials the,first step involved computing 
the equilibrium structure. The starting structure was found from a simulated annealing simulation 
of the isolated N-methylacetamide using our modified empirical force field (Table 1). In this case 
the starting temperature of 300 K was slowly decreased to 5 K during the course of a 15 ps dynam- 
ics run with a time step of 15 fs. The resulting low energy structure ( trans — NMA CC ) was further 
minimized using a steepest descents method until the norm of the forces were < 0.0001. Other 
conformers of NMA were found by an adiabatic mapping of the rotation of the CO and NH methyl 
groups. In the previous ab initio calculations of Mirkin and Krimm, the lowest energy structure for 
fraos-N-Methylacetamide was found to be in a cis-cis configuration (trans — N M A cc ). 16,22 Here 
cis-cis refers to the orientation of the two methyl groups H atom in relation to the CO and NH 
groups. Previous ab initio calculations of Mirkin and Krimm found a 0.01 kcal/mol energy differ- 
ence between trans — NMA CC and trans — N M A tc . 16,22 By adiabatically mapping the rotation of 
the CO methyl group and minimizing using the improved empirical force field, we find the barrier 
of rotation from cis-trans to cic-cis to be 0.23 kcal/mol. In a similar manner, the barrier of rotation 
from cis-cis to trans-cis was found to be 0.65 kcal/mol. The energy difference between the cis-cis 
and trans-cis structures is 0.18 kcal/mol. However, in the current ab initio calculations, we find 
that the lowest energy structure for trans — NMA is in the cis-trans configuration (rrans-NMA ct ), 
as illustrated in figure 1. The cis-cis configuration ( trans — NMA CC ; ) was not a stable minimum 
in the current ab initio calculations and thus it was disregarded for the vibrational analysis. The 
calculated energy difference on the improved empirical surface between cis and trans NMA is 4.5 
kcal/mol in favor of the trans conformer. Ataka and co-workers measured an enthalpic difference 
between the trans and cis conformers to be 1.3 kcal/mol, 13 and previous ab initio results suggest 
an energy difference of 2.5 kcal/mol between these two structures. 16 Our current ab initio calcu- 
lations suggest an energy difference of 2.48 kcal/mol, in complete agreement with the previous 
ab initio calculations of Mirkin and Krimm. 16,15,22 One reason why the empirical Amber-like force 
field could overestimate the energy difference between cis and frans-NMA is a lack of accuracy in 
the description of intramolecular hydrogen bonding present in the cis conformer. 

Due to the length of time required to calculate the ab initio VSCF potential energy surfaces, 
we only calculated the anharmonic vibrational frequencies for the trans — NMA ct configuration 
(Figure 1). The original empirical CC-VSCF, adjusted empirical CC-VSCF, ab initio CC-VSCF, 
and the observed matrix isolated experimental frequencies are compared in Table 2, while the 
harmonic frequencies for all calculations are given in Table 3. The frequency assignments are also 
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given in Tables 2 and 3. The experimental frequencies were obtained in a low temperature nitrogen 
matrix (20K) as reported by Ataka. 13 It is important to mention that in the present calculation, the 
ab initio harmonic and CC-VSCV anharmonic frequencies are not scaled. Indeed, the advantage of 
the CC-VSCF method is that the anharmonic corrections are calculated from the anharmonic part 
of the potential, which is their true physical origin, rather than by a “scaling” procedure, which is 
merely an empirical method based on experience from previous systems. We find the difference in 
the frequencies between the harmonic and CC-VSCF calculations to be as high as 300 cm 1 (Table 
3). One can expect that, in general, an error in the ab initio CC-VSCF method will be on the order 
of 20 cm -1 , based on previous ab initio CC-VSCF calculations, 44-46 and a Density Functional 
calculation. 58 Thus the anharmonicities of the vibrational modes of N-Methylacetamide are quite 
large. 

From Table 2 it is clear that the ab initio CC-VSCF method is far superior to both the original 
and adjusted empirical CC-VSCF calculations. This is particularly true in the case of the amide I, 
II and III modes (1700 to 1265 cm -1 ), which are thought to play a role in intramolecular energy 
redistribution pathways. The improved empirical CC-VSCF calculated frequencies for the amide 
I, II and HI are 1662, 1569 and 1 183 cm -1 respectively. This is compared to the calculated ab 
initio CC-VSCF frequencies of 1751, 1547 and 1283 cm -1 and the experimental frequencies of 
1708, 151 1 and 1265 cm -1 . Both empirical force fields were originally derived for solvated pep- 
tide groups. In solution, the observed N-methylacetamide amide I, II and III frequencies shift to 
1622, 1580 and 1315 cm -1 . 16 The empirical force fields studied in the present work have diffi- 
culty in reproducing matrix isolated NMA vibrational frequencies, particularly in this region. This 
represents a general shortcoming for most empirical force fields used to study high resolution vi- 
brational spectroscopy, i.e. vibrational spectroscopy at a level where the anharmonic effects are 
significant. The fact that experimental data comes from spectroscopy in rare-gas matrices, intro- 
duces an element of uncertainty, since the matrix shift effects can be substantial. However, in a 
recent paper on glycine 46 spectroscopic calculations for ab initio and empirical potentials were 
compared with matrix experiments AND with experiments in superfluid He droplets (where the 
matrix effects are negligible). The superiority of the ab initio potential over the empirical potential 
was found to be even greater when the frequencies were compared with those of the superfluid He 
droplet experiment. 

It is also clear from Table 3 that some vibrational modes of NMA are quite anharmonic, as 
seen by a comparison between the ab initio harmonic and ab initio CC-VSCF frequencies. This is 
particularly true for the stretching frequencies of the NH and methyl groups (3495 to 2900 cm 4 ). 
Because neither the force constants nor frequencies are scaled in the present ab initio calculation, 
the CC-VSCF method has difficulty in reproducing some of the asymmetric and symmetric bend- 
ing modes of the methyl groups (1529 to 1394 cm -1 ). One possible way to improve the accuracy 
of the ab initio calculations is to use a triple-zeta + polarization basis set (TZP). However the ab 
initio CC-VSCF calculation does quite accurately captures the lower frequency modes correspond- 
ing to the in plane and out-of-plane bends and the methyl distortions (< 900 cm 4 ). In this region, 
the empirical expansion method (eq 2. A. 2) becomes less reliable. While the adjusted empirical 
frequencies compare well with both the experimental and anharmonic ab initio frequencies for the 
CHj, asymmetric and symmetric stretches (3000 to 2900 cm,- 1 ), the force field is unable to accu- 
rately describe the the amide I, amide II and amide III frequencies (Table 2). This is however, the 
main strength of the current ab initio CC-VSCF method, the ability to quite accurately describe 
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the anharmonitices of the amide frequencies of a peptide, without resorting to scaling. 

In order to investigate further the degree of anharmonicity in the amide modes, we have defined 
a scalar quantity, the coupling strength (CS), which will indicate the degree of mode-to-mode- 
coupling. Here the coupling strength (CS) for the 1 «- 0 transition, is defined as: 

cs -> h(ij - „°) 


3.1 


In equation 3.1 we divide by the frequency difference between the modes, since in cases of near- 
degeneracy, the effects of a given magnitude of coupling element are much greater. and 'P l - 
represent the ground and first excited state wavefunctions for mode i. The wavefunctions are 
determined self-consistently using equation 2.3, In equation 3.1, is the coupling potential, 
defined by equation 2.B.2 ( ab initio ) and the off-diagonal terms of equation 2.A.2 (empirical). 
Tables 4 and 5 give the matrix elements for the coupling strength C S ^ 0 for both the adjusted 
empirical (Table 4) and ab initio (Table 5) coupling potentials. One would read a symmetric 
matrix element ( CSlj~° ) as representing the coupling strength between modes i and j. 

The largest coupling strength (C5 1<_0 =4.1) calculated using the ab initio coupling potential 
is between the NCH 3 asymmetric stretch and CH$ asymmetric stretch (matrix element 2,5 Table 
5). In the case of the empirical coupling potential, the degree of coupling between the j\CHz 
asymmetric stretch and CH 3 asymmetric stretch is greatly reduced, even though the frequencies 
are nearly degenerate. This is illustrated in figure 2, which is a comparison between the mode- 
to-mode coupling potentials for both the ab initio and the improved empirical force fields. Each 
line in figure 2 represents the coupling potential at a different Qj value for the C H$ asymmetric 
stretching mode at a given value of Q, for the NCH3 asymmetric stretching mode. All normal 
mode coordinates ( Qi ) are in a dimensionless variable as given by equation 2.A.3. It can be seen 
from figure 2 that the ab initio coupling potential is significantly larger than the empirical potential, 
and the ab initio frequencies of these two modes are almost degenerate (14 cm -1 difference). Thus 
the probability of energy transfer between the NCH3 asymmetric stretch and CH3 asymmetric 
stretch is quite large in the ab initio case, 4. 1 , vs 0.0 calculated from the adjusted empirical potential 
energy surface. 

The most interesting finding of the present work is the fact that the amide I mode, at 1708 cm 1 , 
does not couple to other internal modes, as seen by the 8 th row in both coupling strength matrices 
(Tables 4 and 5). This is true for both the empirical and the ab initio coupling potentials. We believe 
that this is a significant finding. In the previous femtosecond nonlinear-infrared spectroscopic 
studies of Hamm and co-workers, the amide I mode of NMAD (deuterated NMA) was found to 
have a vibrational relaxation of 450 fs. 4 This is quite fast, and therefore it was postulated that the 
vibrational energy transfer was internal via an intramolecular energy redistribution (IVR) process. 
Transfer to the solvent should take more time and is expected to be on the order of 10-100 ps. 
Because Hamm and co-workers also saw similar relaxation rates for three small globular proteins, 
they concluded that this relaxation process is essentially the same for the amide I mode of all 
systems,' and is an intrinsic property of the peptide group itself. They proposed that the possible 
candidate modes for energy transfer in NMAD could be the amide IV band at 628 cm 4 , the amide 
III mode at 965 cm -1 or the CO out-of-plane bend at 1044 cm~ l . We had originally thought 
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that the relaxation may be a 2 <— 0 transfer to the amide IH or possibly the amide IV band, as 
there is no appreciable intramolecular coupling for an amide 1 1 <— 0 transitions. The most likely 
reason for this is that the amide I frequency is not near-degenerate to any other intramolecular 
frequencies, and furthermore, the mode-to-mode coupling of amide I is not as strong as other 
calculated intramolecular couplings. Because the coupling strength is weighted by the frequency 
difference, this makes the probability of intramolecular vibrational energy transfer (from 1 <— 0) 
unlikely for the amide I mode in NMA. However, we also calculate the 2 <— 0 coupling element 
between that amide I and other intramolecular modes, and did not find any appreciable coupling 
to either the amide III or the amide IV. Thus, it is possible that the amide I relaxation is due to 
energy transfer to the solvent, or to energy transfer between modes which are not simply expressed 
in a pair-wise fashion. We are currently investigating this further. It therefore still remains unclear 
which modes the experimental vibrational energy of the this mode is transfered to in such a short 
amount of time (450 fs). 

In the case of the amide II mode, however, there is the probability of vibrational energy transfer 
to the other internal modes. For the ab initio case there are two likely intramolecular modes to 
receive energy, the NCH 3 asymmetric bend and the CH 3 asymmetric bend (modes 10 and 1 1). 
The empirical force field also predicts that the amide II couples to the NCH 3 asymmetric bend 
and the CH 3 asymmetric bend, although with different probabilities. However, in the case of the 
empirical force field, we also find significant coupling between the amide II and NCH 3 symmetric 
bend-NC stretch (mode 14), CO out-of-plane bend, CCH 3 rocking (mode 19) and the CO out-of- 
plane bend (mode 21). Figures 3 and 4 illustrate the coupling potentials between the amide II and 
the CH 3 asymmetric bend and the NCH 3 symmetric bend-NC stretch modes (mode 9 coupling to 
modes 1 1 and 14 respectively). In figure 3 it is clear that both force fields couple amide II with the 
CH 3 asymmetric bend. Although the calculated empirical mode-coupled potential energy curves 
are quite a bit larger than the corresponding ab initio calculated mode-coupled potential energy 
curves, the coupling strength is larger in the ab initio case. This is most likely due to the near- 
degeneracy (energy separation of 6 cm -1 ) of the modes compared with an energy separation of 
104 cm -1 in the empirical calculation. 

In figure 4 it is quite clear that the empirical coupling potential greatly overestimates the cou- 
pling between the amide II and the NCH 3 symmetric bend-NC stretch (modes 9 and 14). Similar 
coupling overestimations exist between the other two modes (9 and 19; 9 and 21) as well. In part, 
th strong coupling potentials to these three internal modes may explain why the empirical amide 
II frequency is blue shifted by 22 cm ~ 1 from the ab initio frequency. If one were to simulate an 
intramolecular vibrational energy relaxation process (IVR) between the amide II and other internal 
vibrational modes using an empirical force field, the probability of an erroneous energy transfer 
would exist in this case. In comparison between the empirically calculated and the ab initio cal- 
culated coupling strength, we find that many mode-to-mode couplings are overestimated using the 
empirical force field. It is therefore quite important to compare the ab initio and empirical coupling 
strengths and coupling potentials, in order to accurately gauge possible IVR pathways. 
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IV. Conclusions 


In this paper we address tKe issue of calculating anharmonic vibrational spectroscopy of small 
biological molecules using two different methods to express the coupling between vibrational 
modes. In section I we propose an improved empirical force field, fit to existing spectroscopic data, 
in order to refine the force constants defining the molecule of interest, namely N-methylacetamide 
(NMA). We also calculated the full 30-by-30 pair-wise mode-coupled potential energy surface of 
NMA using the GAMESS electronic structure program with second-order M0ller-Plesset perturba- 
tion theory (MP2). 56 We then calculated the anharmonic vibrational spectrum of N-methylacetamide 
using a Correlation Corrected Vibrational Self-Consistent Field approach, with the effective poten- 
tials defined by either the empirical or the ab initio force fields. Because the resulting CC-VSCF 
ab initio frequencies do not need to be scaled, they can be compared directly to both the empirical 
CC-VSCF and the experimental frequencies (Table 2). 

We find that overall the comparison between the empirical, improved empirical, ab initio and 
experimental frequencies is quite good. In general, however, the ab initio CC-VSCF method il- 
lustrates a significantly better agreement with the experimental data. This particularly true for the 
important and geometrically characteristic modes of amide I, II and III which are found at 1708, 
1511 and 1265 cm -1 . In fact, most of the ab initio CC-VSCF vibrational modes are within 30 
cm -1 of the experimental data, as compared to the improved empirical frequencies which may 
be off by as much as 80 cm~ l . A further comparison between the ab initio harmonic and full 
CC-VSCF calculated frequencies (Table 3) illustrates that most of the ab initio harmonic frequen- 
cies are significantly shifted from their anharmonic counterparts. Because the shift in frequency is 
much larger than the error of the ab initio calculation, this emphasizes the need to utilize a fully 
anharmonic, mode coupled method such as a CC-VSCF calculation when comparing to experi- 
mental vibrational data. All in all, these results suggest the use of ab initio rather than empirical 
potential functions, whenever feasible. 

We further defined a scalar quantity, the coupling strength, in order to gauge the probability of 
vibrational energy transfer processes (1 <— 0 transition). By calculating the intramolecular cou- 
pling strength for both the ab initio and empirical force fields, we find a few surprising results. 
First, the amide I mode does not appear to couple to other internal vibrational modes via a 1 «— 0 
transition. Therefore, it is possible that the experimental results of Hamm and co-workers, 4 which 
suggest that an intrinsic property of the amide I mode is to very quickly (450 fs) transfer vibra- 
tional energy to other internal modes, must come from either coupling directly to the solvent or 
coupling to transitions of lower frequency modes. We also find that the coupling of the amide II 
mode to other internal modes is overestimated by the empirical force field, even if the force field 
is adjusted to fit spectroscopic data. This overestimation could in fact lead to the calculation of 
intramolecular energy redistribution (IVR) pathways for amide II which are incorrect. In general, 
we find that the empirical coupling strength is overestimated for many other modes ranging from 
amide II to amide V. Therefore, a comparison between the coupling strength calculated using an 
empirical and ab initio force fields is a first step in not only benchmarking empirical potentials, 
but also in later developing of biomolecular force fields. Based on this work, we suggest that for 
small biologically relevant molecules it may be better to utilize an ab initio derived potential en- 
ergy surface for spectroscopic applications, as well as for quantitative treatment of intramolecular 
vibrational energy flow. 
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Table I: Molecular Mechanical Parameters 


bond K~ r kcaif(molA 2 ) r e(? X 


Ca-H 

340.0 

1.09 

C-Ca 

317.0 

mESM 

C-O 

570.0 

HUH 

C-N 

490.0 

1.335 

N-H 

434.0 

1.01 

N-Ca 

337.0 

1.449 

Ca-H 

340.0 

1.09 



Angle Parameters 

angle 

Kg kcal /(mol radian 2 ) 9 eg deg 

tt TT — 1 


CO-Ca-H 

40.0 T 

■a 

CO-NH-Ca 

50.0 

329 

Ca-CO-OC 

80.0 

120.4 

Ca-NH-HN 

30.0 

118.04 

Ca-CO-NH 

70.0 

116.60 

H-Ca-H 

38.0 

109.5 

OC-CO-NH 

80.0 

122.9 

CO-NH-HN 

20.0 

120.0 

CO-Ca-H 

40.0 

109.5 

CO-NH-Ca 

50.0 

121.90 

Ca-CO-OC 


120.4 

H-Ca-NH 

40.0 

109.5 

H-Ca-H 

38.0 

109.5 

H-Ca-H 

38.0 

109.5 
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Table I continued: Torsional Parameters 





7 deg 


0.0 

0.0 

0.0 

0.0 

180 

180 

180 

0.0 

90 

0.0 

0.0 


Improper Torsional 


torsion 

K# kcal/(mol radians*) 

<t>e Q deg 

CO-Ca-OC-NH 

100.0 

0.0 

NH-CO-HN-Ca 

45.0 

0.0 














Table 2 : Summary of CC-VSCF Frequencies for trans N-Methyl acetamide 


Empirical 

Adjusted Empirical 

Ab Initio 

Obs. (Matrix) 13 

assignment 

3309 

3310 

3523 

3498 

NH stretch 

2986 

2985 

2993 

2978 

NCH 3 as 

2985 | 

2986 

2985 

3008 

CCH 3 as 

2985 

2986 

3014 

3008 

CH 3 as 

2984 

2985 

2979 

2973 

CH 3 as 

2872 

2873 

2939 

2958 

NCH 3 ss 

2868 

2869 

2950 

2915 

CCH 3 ss 

1676 

1662 

1751 

1708 

amide I 

1598 

1569 

1547 

1511 

amide II 

1466 

1465 

1566 

1472 

NCH 3 asy bend (ab) 

1459 

1459 

1541 

1446 

CH 3 ab 

1461 

1461 

1557 

1432 

CCH 3 ab 

1453 

1461 

1515 

1446 

NCH 3 ab 

1459 

1454 

1468 

1419 

NCH 3 sb NC s 

1388 

1371 

1421 

1370 

CCH 3 ab 

1231 

1183 

1283 

1265 

amide III 

1143 

1148 

1214 

1181 

NH opb, d7 3 rCNs 

1042 

1033 

1184 

1089 

CN s CH 3 t 

1029 

1031 

1119 

1037 

CO opb CCH 3 t 

990 

983 

1083 


CO opb CCH 3 r 

956 

954 

1022 

990 

CH 3 t CC s CO opb 

947 

950 

891 

857 

amide V C H 3 t 

799 

789 

643 

658 

CH 3 r CO ipb 

687 

686 

637 

626 

CCs CO opb 

577 

571 

439 

439 

amide IV C H 3 t 

440 

402 

487 

391 

CO opb CCH 3 t 

302 

281 

299 

279 

CNC d CCN d 

227 

234 

245 


CNC d CCN d CCH 3 d 

372 

381 

266 


NCH 3 d 

398 

398 

166 


CCH 3 d 


s=stretch, as=asymeteric stretch, ss=symmetric stretch, d=distortion, crocking, ab=asymmetric bend, 


sb=symmetric bend, ipb= in-plane bend, opb=out-of-plane bend 
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Table 3 : Summary of Harmonic Frequencies for frans N-Methylacetamide 


Empirical 

Adjusted Empirical 

Ab Initio 

Obs. (Matrix) 13 

assign 

3316 

3316 

3751 

3498 

NH stretch 

2992 

2992 

3261 

2978 

NCH 3 as 

2992 

2992 

3248 

3008 

CCH 3 as 

2992 

2992 

3245 

3008 

CH 3 as 

2991 

2991 

3225 

2973 

CH 3 as 

2878 

2878 

3133 

2958 

NCH 3 ss 

2875 

2875 

3124 

2915 

CCH 3 ss 

1679 

1666 

1780 

1708 

amide I 

1597 

1567 

1584 

1511 

amide II 

1466 

1464 

1548 

1472 

NCH 3 asy bend (ab) 

1457 

1457 

1526 

1446 

CH 3 ab 

1453 

1453 

1520 

1432 

CCH 3 ab 

1453 

1453 

1508 

1446 

NCH 3 ab 

1447 

1447 

1484 

1419 

NCH 3 sb NC s 

1380 

1359 

1438 

1370 

CCH 3 ab 

1218 

1169 

1309 

1265 

amide HI 

1134 

1132 

1208 

1181 

NH opb, CH 3 r CNs 

1021 

1010 

1173 

1089 

CN s CH 3 t 

984 

984 

1133 

1037 

CO opb CCH 3 t 

968 

952 

1074 


CO opb CCH 3 t 

927 

922 

1020 

990 

CH 3 t CC s CO opb 

903 

903 

890 

857 

amide V CH 3 r 

793 

780 

634 

658 

CH 3 t CO ipb 

661 

659 

624 

626 

CCs CO opb 

575 

570 

428 

439 

amide IV CH 3 r 

438 

399 

362 

391 

CO opb CCH 3 t 

295 

272 

270 

279 

CNC d CCN d 

217 

219 

153 


CNC d CCN d CCH 3 d 

351 

351 

53 


NCH 3 d 

361 

360 

166 


CCH 3 d 

— " « » i . • i. _ 


s=stretch, as=asymeteric stretch, ss=symmetric stretch, d=distortion, r=roc king, ab— asymmetric bend, 
sb=symmetric bend, ipb= in-plane bend, opb=out-of-plane bend 
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CNC Angle 120.7(121.95) 
NCC Angle 115(117) 


CNCO Torsion -3.62 (0.08) 
HNCO Torsion -173.34 (—179.9) 


Figure 1 : The optimized equilibrium structure of fra/is-N-Methylacetamide. Bond lengths (in A) 
and angles (in degrees) are given for the ab initio MP2/DZP and empirical (parenthesis) optimized 
structure of NMA. 



Normal Mode Coordinates, NCH3 as 



Figure 2: The coupling potential between the NCH 3 asymmetric stretch and the CH 3 asymmetric 
stretch (modes 2 and 5) for the ab initio (top panel) and the empirical (bottom panel) force fields. 
The normal mode coordinates are in a dimensionless variable as given by equation 2.A.3. The 
coupling potential is given in Hartrees. For comparison the potentials are given on the same scale. 






Figure 3: The coupling potential between the amide II and the CH 3 ab (modes 9 and 1 1 ) for the ab 
initio (top panel) and the empirical (bottom panel) force fields. The normal mode coordinates are 
ip a dimensionless variable as given by equation 2. A. 3. The coupling potential is given in Hartrees. 
For comparison the potentials are given on the same scale. 
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Figure 4: The coupling coupling potential between the Amide II and the NCH 3 sb NC s (modes 9 
and 14) for the ab initio (top panel) and the empirical (bottom panel) force fields. The normal mode 
coordinates are in a dimensionless variable as given by equation 2.A.3. The coupling potential is 
given in Hartrees. 



Table 4: Coupling Strength for the Adjusted Empirical Force Field 



Table 5: Coupling Strength for the Ab Initio Force Field 
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;=stretch, as=asymeteric stretch, ss=symmetric stretch, distortion, r==rocking,ab=asyrnmetric bend, sb=symmetrtc- 


